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This paper provides insight on the economic inefficiency of the classical merit-order dispatch in electricity 
markets with uncertain supply. For this, we consider a power system whose operation is driven by a two-stage 
electricity market, with a forward and a real-time market. We analyze two different clearing mechanisms: 
a conventional one, whereby the forward and the balancing markets are independently cleared following a 
merit order, and a stochastic one, whereby both market stages are co-optimized with a view to minimizing 
the expected aggregate system operating cost. We first derive analytical formulae to determine the dispatch 
rule prompted by the co-optimized two-stage market for a stylized power system with flexible, inflexible and 
stochastic power generation and infinite transmission capacity. This exercise sheds light on the conditions for 
the stochastic market-clearing mechanism to break the merit order. We then introduce and characterize two 
enhanced variants of the conventional two-stage market that result in either price-consistent or cost-efficient 
merit-order dispatch solutions, respectively. The first of these variants corresponds to a conventional two- 
stage market that allows for virtual bidding, while the second requires that the stochastic power production 
be centrally dispatched. Finally, we discuss the practical implications of our analytical results and illustrate 
our conclusions through examples. 
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1. Introduction 

Electricity markets are typically arranged as sequences of exchanges or pools, where producers 
and possibly consumers submit offers and bids specifying the amount of electricity they are willing 
to deliver to or withdraw from the network and at what unit price. During the market-clearing 
process, a market operator determines the optimal dispatch by accepting a subset of the submitted 
offers for electricity production and, possibly, bids for consumption. Furthermore, a market price 
or a set of prices, one for each node across the electricity network, is set so that the dispatched 
generation and consumption blocks are profitable according to the respective offer and bid prices. 

Generally, electricity markets comprise different floors for trading electricity arranged in a sequen¬ 
tial fashion up to the time when electricity is delivered. The presence of a real-time or balancing 
market is necessary as electricity is a non-storable commodity, and since imbalances between sup¬ 
ply and offtake result in deviations of the system frequency that can harm machines and appliances 
connected to the network. On the other hand, market floors that clear hours ahead of electricity 
delivery are also needed to guarantee the participation of units with slower response time, e.g., 
nuclear power plants. Although market structures and regulations vary significantly from country 
to country, a common trait is the presence of at least a day-ahead market stage (besides a real¬ 
time one) clearing around noon on the day prior to the delivery of electricity. In m any Europea n 


countries, day-ahead markets account for the bulk of the total trading of electricity (j Weber 


20101 


Although other market stages may exist, we consider in the following a two-stage market compris¬ 
ing a forward (day-ahead) and a real-time floor as an abstraction of electricity markets comprising 
different market stages. 

Traditionally, electricity markets are cleared in a sequential manner and independently from 
one another, i.e., without considering the impact of the market-clearing decision on the operation 
of future markets stages. In such a conventional arrangement, each market stage is cleared by 
dispatching producers and consumers in order according to increasing marginal costs and decreasing 


marginal utilities, respectively. Under the assnmption that producers and consnmers bid their 
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“true” marginal costs and benefits, such a dispatch of the forward market based on the merit-order 
principle results in schedules that maximize the social welfare exclusively for that market stage. 
However, there is no guarantee that it maximizes the total social welfare, i.e., aggregated across 
sequential market stages, in the long run. Indeed, the application of the merit-order principle may 
leave fast-ramping units out of the forward dispatch and potentially result in a lack of flexibility, 
and hence inefficiency, at the real-time stage. 

Argnably, the suboptimality of a merit-order based forward dispatch is exacerbated by the 
increasing penetration of partly-predictable renewables such as wind and solar. These generation 
sources are characterized by a zero (or near-zero) production cost per nnit and their dispatch is 
often prioritized by market operators. As a result, owners of renewable production facilities submit 
price-inelastic forward offers, where a certaiii quant ity, typically a point forecast of future gener¬ 


ation, is offered at zero price (IMorales et al 


2ni4al l. In a merit-order based forward dispatch, a 


large penetration of renewables may, besides pushing flexible resources out of the market, increase 
the need for flexibility at the balancing stage, which resnlts from their uncertain nature. 

The attention received by the impact of uncertain renewable generation on electricity markets 
has grown in the technical li terature along with t heir actual deployment in power systems. Numer¬ 


ous papers, see for example 


Jonssonetal 


(120101 ). have assessed the downward pressnre exercised 


by uncertain renewable power on electricity market prices, i.e., the so-called merit-order effect. 
Among the most notable consequences of the merit-order effect is the emergence of negative prices, 
which are caused by the combination of large volnmes of zero-price offers from renewable suppliers 


and t 
units. 


lermal producer s that are willing to incur occasional losses to avoid wear-and-tear of their 


Hildmann et al. 


(j2014 1 claims that the removal of feed-in-tariff schemes wonld partially solve 


this problem by inciting owners of renewable generation facilities to internalize forecasting-error 
(balancing) costs in their offer, hence resulting in positive marginal costs of renewables. However, 
balancing costs are hard to estimate and vary from hour to hour. 

Other works have foc used on centralized solutions for the coordin ated clearing of the forward and 


real-time market stages (jPritchard et al. 


2010 


Morales et al 


20121 ). In these proposals, the forward 
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dispatch is determined by a market operator that makes use of two-stage stochastic programming 
to account for the balancing costs due to the uncertain supply and minimize the expected aggregate 


system cost. While stochastic dispatch models do not comply with the merit-order at t 


le forward 


stage, they guarantee both reven ue adequacy for the market operator (|Pritchard et al. 


cost-recovery for the producers (IMorales et al. 


3 


201011 and 


20121 1 in expectation. Furthermore, they improve 


the long-run social welfare. An alternative solution where the merit-order at the for ward market 
is only broken for providers of uncertain supply is proposed in 


Morales et al. 


(j2014bl l. That work 


shows that part of the improvement in expected social welfare achievable by the stochastic dispatch 
can be captured by a classical forward market where the dispatch of uncertain renewable generators 
is forced by the market operator. In general, though, the problem of determining the optimal 
forward offer or dispatch for stochastic power suppliers has no c losed form so l ution in terms of a 
simple statistic of the forecast distribution of uncertain supply. 


Zavala et al 


(|2015fl investigates 


stochastic dispatch schemes with and £2 penalization of imbalances where the optimal dispatch 
converges to the conditional median and mean of the uncertain supply. However, a generalization of 
these properties is not possible, as the underlying penalty assumptions may clash with the offering 
preferences of the flexible producers, which ultimately drive the penalization of the imbalances at 
the real-time market. 

This paper aims to shed light on the impact that preserving the merit order in forward electricity 
markets with uncertain supply may have on market efficiency. By market efficiency we mean the 
ability of the market to minimize the expected aggregate system operating cost (cost efficiency) 
and to deliver a set of prices such that the forward price equals the expectation of the real-time 


price (price cons istency) — a number of aut 


riors 


market settings (jBessembinder and Lemmon 


rave underlined t 


2002 


Kave et al 


19901. 


e benefits of 

price-consistent 

Zavala et al. 

2015 

). For this 


purpose, we construct mathematical models for four different types of two-stage markets, namely; 
1. A market that is price-consistent and guarantees maximum cost-efficiency. A market with such 
properties is achieved by co-optimizing the forward dispatch and the real-time re-dispatch 
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through the use of stochastic programming. We show that this market produces, however, dis¬ 
patch solutions that break the merit order. Furthermore, the practical implementation of this 
market is not without its challenges and problems in terms of revenue adequacy, cost recovery. 


arbitrariness in the probabilistic characterization of the uncertain supply, etc. (IMorales et al. 


2014b|). Therefore, we just use it here as an “ideal” benchmark that theoretically achieves the 


highest market efficiency. 

2. A market that follows the merit order, but that is, in general, price inconsistent and cost- 
inefhcient. This is the case of a conventional two-stage market, where the forward and the 
real-time settlements are not co-optimized, and where the uncertain supply is systematically 
dispatched to a certain statistic of its forecast probability distribution (most commonly, the 
conditional expectation). 

3. A price-consistent market that preserves the merit order. We construct this market from the 
conventional two-stage market described in point 2 above, by introducing a risk-neutral virtual 
bidder that arbitrages between the forward and the real-time markets. We show that, in order 
to ensure price consistency, such a market may have to give up on cost efficiency. 

4. A market that renders the most cost-efficient dispatch among those that respect the merit 
order. The practical translation of this market is that of a conventional two-stage market, like 
the one described in point 2 above, in which the uncertain supply is centrally dispatched by a 
non-profit, all-knowing organization such as an Independent System Operator. We show that, 
in order for this market to ensure maximum cost-efficiency, while complying with the merit 
order, it may have to give up on price consistency. 


Unlike other works th a t rely on comput a tional simulat i on for the ir analysis, e.g.. 


Bouffard and Galianal (j2008f) . 


Khazaei et al 


mm . 


Morales et al 


(201^ 


2014bD, we derive closed- 


form solutions to the mathematical models describing these four types of markets for a stylized 
power system with flexible, inflexible and stochastic power generation and infinite transmission 
capacity. This exercise allows us to characterize the dispatch solutions prompted by these markets 
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and identify conditions for their equivalence or dissimilarity. We accompany this analytical insight 
with a meaningful discussion on the practical implications of our results and illustrate our main 
conclusions through examples. 

The structure of the paper is the following. Sections [2] and [3] deal with markets 1 and 2, respec¬ 
tively. More specifically, we define and formulate the conventional and the stochastic two-stage 
electricity markets for a stylized power system and provide the closed forms of the dispatch rules 
that each of these markets induce. In these two sections we also introduce some important concepts 
that are repeatedly used throughout the paper. In Section 01 we focus on market 3 and provide 
conditions under which the conventional and the stochastic market-clearing models are equivalent 
in the case that virtual bidding is allowed. Section [5] deals with market 4, that is, with the case of 
a conventional two-stage market in which the stochastic power production is centrally dispatched 
with the aim of minimizing the expected aggregate system operating cost. The study of this mar¬ 
ket setting allows us to identify conditions under which a conventional two-stage market, even if 
price-consistent, does not deliver the most cost-efficient merit-order dispatch. Finally, conclusions 
are drawn in Sectional 

2. Conventional or InefFicient Two-stage Market (ConvM) 

We first formulate the model of the conventional two-stage market, where the operation of the 
forward and the balancing markets is not co-optimized. Each market, therefore, attempts to min¬ 
imize operating costs independently. We describe such a two-stage market as inefficient, because 
it results in higher operating costs in the long run. 

The market models that we introduce throughout the paper are all tailored to the stylized power 
system described below. 

Definition 1 (Stylized power system). Our stylized power system has infinite transmission 
capacity and consists of inflexible, flexible and stochastic power generation technologies with capac¬ 
ities Pj > 0, Pp > 0, and pyp > 0, in that order. We denote the marginal cost of the flexible and 
inflexible generating capacity by > 0 and C/ > 0, respectively, and assume that the marginal 
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cost of the stochastic power production is zero. Furthermore, we use to represent the incre¬ 
mental cost incurred by the flexible generating capacity for marginally increasing its production 
for balancing (upward regulation) and to denote the incremental utility that it obtains from 
marginally decreasing its production for balancing (downward regulation). The demand I is inelas¬ 
tic and known with certainty with a cost of involuntary curtailment denoted by v. In this stylized 
power system, it holds that v > Cp > Cp > Cp > 0. Finally, the power production from renewable 
sources is characterized as a random variable W defined on some probability space Let 

F(-) and /(•) denote the cumulative distribution function and the probability density function of 
IF, respectively. To comfortably deal with the case F{0) >0 in our mathematical derivations, we 
use the generalized inverse cumulative distribution function F~^{a) = inf{x € M : F(x) > a}. 

The aim of the forward market is to determine the dispatch of flexible power producers (pf), 
inflexible power producers (pi), and stochastic power producers (pw) that minimizes (forward) 
system operating costs, that is, 


Minimize Cjpj + CfPf 

PF,PI,PW 

(la) 

s.t. Pf +Pi+Pw - 1 = 0 ■ 

(lb) 

0<Pf<Pf, 

(Ic) 

0 < Pi < Pi, 

(Id) 

0 < Pw < Pw 

(le) 


Equation ffTbl) enforces the power balance. We denote the Lagrange multiplier associated with 
the power balance equation by Xf which defines the forward electricity price. The set of inequal¬ 
ities (fTc)l - (fTel) imposes upper and lower bounds on the dispatch of the different power producers. 
We explicitly indicate parameter p^f, which stands for the power production that is expected from 
the stochastic power producers, in (fTe)l . We do so to note that, typically, the amount of stochas- 
tic power production that can be cleared in the forward market is capped to this expec tation 


flBouffard and Galiana 


2008, 


Cadre and Didier 


2014 


Qggioni et al 


2014 


Zavala et al. 


201 ,^ . 


The following proposition provides the optimal solution to problem m- 
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Proposition 1 (The merit-order dispatch solution). Consider the stylized power system 
described in Definition {![ where, in addition, it holds that Cj < Cp- Optimization problem © 
prompts the following dispatch rule: 


Rule # 

Pw Pi 

Pf 

applies if 

1. 

1 0 

0 

0 < / < pw 


2. 

1 

o 

Pw < 1 < Pw + Pi 

3. 

Pw Pi 

1 

1 

Pw +Pi <l<pw +Pi+Pf 



p .7 , 




4- 

27 


pw + Pj - ■ 




Gomez-Exposito et al 


(l2008l. Chapter 5)) and therefore, its 


This result is well known (see, e.g., 
proof is omitted here. 

Note that values of load I > pw +Pi+Pf render problem ([1]) infeasible, because we have not 
considered the possibility of shedding load in the forward market. 

The conventional market produces a forward dispatch whereby power production is cleared fol¬ 
lowing the so-called merit order, i.e., the power plants with the lowest marginal costs are dispatched 
first, until the system demand is satished. In this case, the dispatch solution prompted by this 
market follows from the intersection of the (forward) supply cost function of the system with the 
marginal utility demand curve. The supply cost function is built by sorting the marginal cost func¬ 
tions of the individual power plants in increasing order. Consequently, all the production with a 
marginal cost lower than the one determined by the said intersection is dispatched. 

Optimization problem ([T|) results, therefore, in the forward dispatch {pf,Pi,Pw) given by Propo¬ 
sition [H Since the stochastic power production cannot be perfectly predicted, the (random) power 
imbalance W — pw is to be covered in the balancing market. For this purpose, flexible power plants 
can be re-dispatched and/or the amount of load and stochastic power production can be curtailed. 

Suppose a specific realization W(uj), a; G id, of the stochastic power production. The balancing 
market determines the most economical vector of re-dispatch actions that accommodates the power 
imbalance —pw, that is. 


Minimize vs{uj) + c'fpf{uj) — CpPpi^) 


(2a) 
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s{uj) + pf{uj) -Pf{uj) + Apw{uj) = 0 : X^{u}) , 

(2b) 

0<pf(uj) <Pf , 

(2c) 

0<pf(uj) <Pf-Pf , 

(2d) 

0 ^ Pw T ^Pw (^) ^ hP (w) , 

(2e) 

0 < s(cj) < 1 , 

(2f) 


Equation (l2bl) ensures that the power system is brought to balance by deploying upward or 
downward regulation from the flexible power unit, i.e., pfiuj) or pp(uj), respectively; curtailing 
load s(w) and/or curtailing stochastic power production, which is given by —pw — ^Pw{'^)- 

The Lagrange multiplier A*'(w) defines the marginal price that clears the balancing market. The 
family of inequalities set limits on the amount of downward and upward regulation that 

the flexible power unit can provide, (I2B and ()2dll . respectively; the amount of stochastic power 
production that can be curtailed ([^ . and the amount of load that can be shed ([2^ . 

The total expected cost of operating the power system is given by C/P/ + Cppp +U{pj,pp,p]y), 
where U{pi,pF,Pw) represents the expected balancing cost computed in the proposition below. 

Proposition 2 (Expected balancing cost). Consider the stylized power system described in 
Definition [7] with given forward dispatch quantities pi,pF,Pw- The expected balancing cost 
U{Pi,Pf,Pw) is computed as 

rPF+Pw-PF rPw rPF+Pw 

C^{Pi-,Pf-,Pw) = v / F{oj)duj + c'f / F{ixi)duj + c~f / F{uj)duj — cfpF (3) 

0 '^pf'^pw~pf Pw 

The proof of this proposition is included in Appendix El Note that the expected balancing cost 
is actually independent ofp/. 

3. Stochastic or Efficient Two-stage Market (StoM) 

We now build the model of a two-stage market where the operation of the forward and the balancing 
markets is co-optimized. To this aim, one just needs to replace optimization problem ([T|) with an 
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alternative market-clearing mechanism that seeks to minimize the expected total system operating 
cost, namely: 


Minimize CjPj + CfPf+ / {vs{uj) 

,pw'FtP)’PpP)PpwP)<^P) dn 

-1- c'^pp{uj) - CpPp{u})) f{uj)duj 

(4a) 

Pf+Pi+Pw-}-=0:U 


(4b) 

s{u}) + p+{uj) -pp{u;) + Apvv(w) = 0 : 

Vw G 0 

(4c) 

0 <Pf <Pf 


(4d) 

0 < Pi < Pi 


(4e) 

0<Pp{u})<Pf, Vw € 


(4f) 

0 <pp{uj) <Pf — Pf, Vo; G n 


(4g) 

0 <Pw + Apvc(w) < IF(w) : 7(w),7(w), Vcu G 0 


(4h) 

0 < s(w) <1, Vo; G n 


(4i) 


where the expectation of the balancing cost (^vs(uj) + CpPp(uj) — CpPp{u)) f{uj)duj is taken over 
the probability space on which the stochastic power production W is defined. Note that, 

for ease of notation, we write y(w) instead of y{W{u})). 

Problem ([4]) computes the optimal forward dispatch {pf,Pi,Pw) by taking into account 
the potential cost of the subsequent re-dispatch of the system that is induced by the ran¬ 
dom power imbalance W — pw Ideally, problem (jl]) also provides the optimal re-dispatch rule 
(p^(w),p^(a;), Apvv(w), 5 ( 0 ;)) that guarantees, by enforcing and (|lj)-(l4i]), the power balance 
for any possible outcome W{ijo) of the random variable W. 

We note that, unlike ([I]), the forward market dH does not need to arbitrarily cap the dispatch pw 
of stochastic power production, because in Q the optimization of pw is driven by the probabilistic 
characterization of random variable W, which is naturally bounded and nonnegative. By the same 
token, we do not need to impose that pw > 0, since c'^>CF>Cp> 0. 

Finally, the balancing stage of this market is also modeled by m, but with the optimal forward 
dispatch {pf,Pi,Pw) given by (H]). 

We now define some relevant concepts that will be used in the remaining part of this paper. 



Morales and Pineda: On the Inefficiency of the Merit Order in Forward Markets with Uncertain Supply 
Article submitted to Operations Research-, manuscript no. (Please, provide the manuscript number!) 


11 


Definition 2. A two-stage market such as (HI)-© and (HI)-© is said to be price consistent if the 
forward price is equal to the expected value of the real-time price, that is, 


Xf = En [A'(a;)] = / X\u;)f{u)du 

Jn 


in the case of the conventional two-stage market m-m, and 

U = Eq [n\u))] = [ v^{u})f{uj)di 

Jn 


(5) 


w 


( 6 ) 


in the case of the co-optimized two-stage market (HI)-©- 


Zavala et al 


( 2015 1. Definition [2] allows us to 


We have taken the term price consistency from 
formulate the following proposition. 

Proposition 3 (Price consistency of the stochastic two-stage market). The two-stage 
stochastic market @-( 121 ) is price consistent. 

Proof. The optimality conditions of problem (|3]) imply that the derivative of the Lagrangian 
with respect to pw and Apvr(<^) are equal to 0 at the optimal solution, i.e., 

^ = 0 ^ U + (-7(w) + 7(‘^)) du = 0 

= 0 i^*'(a;)/(a;) — 7(w)-1-7(0;) = 0 ,Vw € D 

“ □ 


i/t = 


dC 


Jn 


(7) 


9Apvi.r(cj) 

where 7 ( 0 ;) and 7 ( 0 ;) are the dual variables of constraints (l4h|) . 


Definition 3. Consider the stochastic market-clearing mechanism ([3]). This mechanism is said to 
break or violate the merit order when a conventional generating unit is dispatched in the forward 
market ahead of a power generating unit with a lower marginal cost. 

Note that this latter definition only applies among conventional power plants and leaves out 
stochastic power generating units. In the particular case of the stylized power system introduced 
in Definition [H the merit order is violated if pi < pj and pp > 0 (with cj < Cp) or if pp < pp and 
Vt > 0 (with cf < Ct). 


Morales et al 


( 120121 ) and 


Morales et al 


(l2014bl) provide examples of situations in which the 


generation dispatch {pp,Pi,Pw) that is solution to the stochastic market-clearing mechanism ([3|) 
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violates the merit order. In this paper, we aim at providing analytical insight on the conditions 
under which this occurs. Our analysis relies on the following theorem, which is indeed one of the 
major results of our work. 

Theorem 1 (The stochastic dispatch rule). Consider the stylized power system described in 
DefinitionUl Now define the constants 



^ 2 :=min I: {v - c'^) F {I - pp) + CpF{l) > Cj (9) 

/ 3 :=min I: Cp + {v - cj,) F {I-pp) - Cp {1 - F{1)) > cp, (10) 

^4 := min I: (v - cf) F {I-pp) + CpF{l) > Cj (11) 

i>ii+Pp 

to which we assign an infinite value in those cases where the corresponding minimization problem 
is infeasible. The efficient two-stage market Q prompts the following dispatch rule, which we refer 
to as the stochastic dispatch rule; 


Rule # 

Pw 

Pi 

Pf 

applies if 

1. 

1 

0 

0 


0<l<l2 



2. 

h 

l-h 

0 


h < 1 '^Pj + h 



3. 

i-Pi 

Pi 

0 

h > h 

Pj + I 2 <l <Pi 5-h 



4- 

h 

Pi 

1 — h —Pi 


Pj + h <l <Pf +Pi 5-h 


5. 

i-pp-Pi 

Pi 

Pf 


l>PF5-Pi + h 



6 . 

1 

0 

0 


0<l<h 

7 . 

h 

0 

l-h 



VI 

VI 


8. 

h 

l-h 

h — h 



13 <l ^Pl -\-13 






I 3 < h 



9. 

h 

Pi 

1 - h -Pi 



h +Pi < 1 <Pi+PF + h 

10. 

I-Pi-Pf 

Pi 

Pf 

h < h 


1 >pj -\-pp-\-h 

11. 

h 

0 

l-h 



h <l <Pf + h 

12. 

1 — Pf 

0 

Pf 



314 

Pf + h <l <h 

13. 

U - Pf 

l-U 

Pf 


I 3 > h +Pf 

[ci < v) 

/d /d ~\~Pj 

U- 

1 

1 

Pi 

Pf 




I>l4+Pj 

15. 

1 — Pf 

0 

Pf 



$l4 

l>Pp3-h 
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The proof of this theorem is given in Appendix [Bj 

Theorem [1] can be used to identify cases in which the stochastic market-clearing mechanism Q 
breaks the merit order. For example, consider a power system whose ability to cope with uncertain 
supply stems from a small amount of flexible power capacity that is comparatively cheap {cp < 
Cl -C v) and that can provide upward regulation at almost no extra cost {cp < Cp). Under these 
conditions, the block of dispatch rules 1-5 {li > I 2 ) in Theorem [1] may apply, resulting in the 
inflexible generation technology being dispatched over the cheaper flexible one. Intuitively, it makes 
economic sense, under these circumstances, to withhold flexible capacity from the forward market 
to have it available for upward regulation in real time at almost no extra cost (bear in mind that 
the alternative would be to curtail load at the very high cost v). 

The mirrored case would be that of a power system with a large amount of flexible, but com¬ 
paratively expensive power capacity (c/ <cp <^v) that is able to provide downward regulation at 
nearly no extra cost {cp < Cp). In this situation, the block of dispatch rules 6-15 {li < I 2 ) would 
hold and the flexible capacity would be dispatched over the cheaper inflexible one. The intuition 
here is that it is profitable to commit flexible capacity in the forward market to have it available 
for accommodating surplus of stochastic power production in real time at nearly no extra cost. 

In the sequel, we focus our analysis to those cases in which cj < Cp. This is a fairly common 
characteristic of many power systems, where, for instance, the base load is mainly supplied by 
nuclear and large coal-fired power plants, while balancing is mostly provided by gas-fired power 
units. Notice that the fact that C/ < Cp precludes the stochastic dispatch rule 15 in Theorem [H 
Next we provide results that are directly derived from Theorem [1] for relevant special cases. We 
state these results in the form of corollaries. 

Corollary I (Free downward regulation). Consider the power system deserihed in Defini- 
tionUl where, in addition, we have that 0 < C/ < Cp, Cp — c~p = (free downward regulation), and 
F'“^(0) = 0. The stoehastic dispatch rule simplifies to: 
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Rule # 

Pw 

Pi 

Pf 

applies if 

7. 

0 

0 

1 


0 < / ^ /s 

8. 

0 

l-h 

h 

VI 

/s < ^ < Pj + ^3 

9. 

0 

Pi 

l-Pi 


I3+P1 < 1 <Pi+Pf 

10. 

l-Pi-P f 

Pi 

Pf 


1 >Pi+Pf 

11. 

0 

0 

1 


0 <l <Pp 

12. 

1 -Pf 

0 

Pf 


Pp 1 I4 

13. 

I4, pp 

l-u 

Pf 

I3 >Pf 

I4 ^ ^ Pj 

14 . 

1 -Pf-Pi 

Pi 

Pf 


l>h+Pi 


Proof. The proof of this corollary directly follows from Theorem [T] and the fact that if cp — cf = 
0, then/i = F“^(0) = 0 </ 2 - D 

The case described in this corollary is, perhaps, the most paradigmatic example of a power 
system in which maximum efficiency is to be achieved by breaking the merit order in the forward 
market. Indeed, if the provision of downward regulation does not entail any extra cost to the 
system (situation that we describe as “free downward regulation”), the more expensive flexible 
power capacity should be always dispatched first, even before the cheaper inffexible one, in the hope 
that the dispatched flexible capacity can be de-committed to accommodate the eventual stochastic 
power production. One can think of the provision of downward regulation as a sort of arbitrage 
whereby the flexible generation technology is, in the e nd, requested not to produce the amount 


of power that it was scheduled to supply in advance (jPritchard et al 


201(111 . Consequently, one 


can argue about the reasons why the cost of downward regulation should be different from zero, 
that is, Cp — cf ^ Q (such as the potential existence of nonconvex costs that are not captured in 
self-commitment based electricity markets). 

We now deal with the mirrored case of free upward regulation. 


Corollary 2 (Free upward regulation). Consider the power system described in DefinitionUi 
where, besides, we have that t) < Ci < Cp and Cp — cf = 0 (free upward regulation). The stochastic 
dispatch rule boils down to: 
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Rule # 

Pw 

PI 

Pf 

applies if 

1 . 

1 

0 

0 

0 ^ 1 ^ I 2 

2 . 


l-h 

0 

^2 < ^ ^ P/ + ^2 

3. 

l-Pi 

Pi 

0 

Pi + h<l<Pi +Pw 

4- 

Pw 

Pi 

^ “ Pw — Pi 

Pi +Pw < 1 <Pf +Pi +Pw 

5. 

1-Pf-Pi 

Pi 

Pf 

1>Pf+Pi+Pw 


Proof. Again the proof of this corollary directly follows from Theorem [U by noting that Cp > C/ 
and Cp — 0 ^ = 0 implies li = F~^{1) =Pyy. Since I 2 < p^, then I 2 <li. □ 

In the case that the provision of upward regulation does not impose any extra cost on the 
system (situation that we refer to as “free upward regulation”), the stochastic market-clearing 
mechanism ([3|) prompts forward dispatch solutions that respect the merit order in the forward 
market (that is, the more expensive flexible power capacity is only dispatched differently from zero 
for levels of load at which the inflexible generation technology is requested to produce at maximum 
capacity). As we will show later, this family of dispatch solutions (namely, those honouring the 
merit order) can be reproduced, at least partially, by enhanced forms of the conventional two-stage 
market. 

Corollaries [Hand [2] define two extreme cases of a broader family of power systems characterized by 
asymmetric costs for balancing power. In a power system where the cost differential for downward 
regulation {cp — cf) is sufficiently lower than the cost differential for upward regulation {cf — Cp), 
the conventional two-stage market will prove to be inefficient, because the most economical way 
of operating the power system is that for which the dispatch of the flexible generating capacity is 
prioritized over the commitment of the cheaper inflexible one. In the opposite case, where Cp — cf 
is enough greater than cf — Cp, the most cost-efficient forward dispatch complies with the merit 
order and, therefore, could still be induced by the conventional two-stage market. Nevertheless, 
one could argue for the case cf — Cp > Cp — cf to be the general rule, insofar as the provision of 
upward regulation involves generating more energy than scheduled, while the provision of downward 
regulation entails not to honour a forward contract. 
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We now complement Corollaries [T] and [2] with the following result, which pertains to the case for 
which the marginal costs of upward and downward regulation are not marked up with respect to 
the forward production costs, that is, Cp — Cp = Cp — Cp = 0. 

Corollary 3 (Free upward and downward regnlation). Consider the power system 
described in Definitions^ where, besides, we have that 0 < C/ < and Cp — = Cp — Cp = 0 (free 

upward and downward regulation). The stochastic dispatch rule reduces to: 


Rule # 

PW + PF 

pi 

applies if 

1. 

1 

0 

0<l<l2 

2. 

h 

l-h 

^2 ^ ^ Pj I 2 

3. 

i-Pi 

Pi 

Pj + I 2 <i 


with pw > 0 and pp > 0. 


Proof. This result can be proved by noting that if — cj = = 0, then h = h = U- D 

Corollary ([3]) describes a power system for which dispatch and re-dispatch actions are equally 
costly. In this case, there is always a forward dispatch solution that satisfies the merit order (by 
just setting = 0 in rules 1 and 2). As we will see later, this makes it possible to analyze enhanced 
variants of the conventional two-stage market that close the cost-efficiency gap with respect to the 
ideal stochastic market-clearing mechanism dH). 

In the stochastic dispatch rule given by Theorem [H the capacities of the flexible and inflexible 
generation technologies play roles that are as important as their marginal generation costs. In this 
vein. Corollary [3] below provides two relevant results. The hrst one pertains to the case of a power 
system where there is enough flexible power capacity to cover any potential lack of stochastic power 
production, that is, pp > pyp. In an abuse of terminology, we will refer to this instance as “capacity 
adequate power system”. The second result constitutes a further simplihcation of the stochastic 
dispatch rule that is possible when, in addition, the range of system loads does not exceed the 
capacity of the inflexible generation technology, that is, l<Pj. 


Corollary 4 (Capacity adequate power system). Consider the power system described in 
Definitions^ where, in addition, we have that cj < Cp and pp >Pw The stochastic dispatch rule 
simplifies to: 
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Rule # 

Pw 

Pi 

Pf 

applies if 

1. 

1 

0 

0 


0<l<l2 

2. 

I 2 

l-h 

0 


I 2 ^ ^ ^ Pi ^2 

3. 

l-Pi 

Pi 

0 

h > I 2 

Pi + l2<l<Pi + h 

1 

h 

Pi 

l — h—Pj 


Pi + h < 1 <Pf +Pi + h 

5. 

l-Pp-Pi 

Pi 

Pf 


^ > Pf 'h Pi + h 

6. 

1 

0 

0 


0<l<h 

7. 

h 

0 

l-h 


CO 

VI 

VI 

8 . 

h 

l-h 

h — h 

h < I 2 

13 <l ^Pl 13 

9. 

h 

Pi 

l-h-P i 


h+Pi <l<Pi+PF + h 

10. 

l-Pi-P f 

Pi 

Pf 


l>Pi+Pp + h 


where: 

= and I 3 = max , F~^ ^1 — ^ ^ with li < <pp. 

If, in addition, we have that I <pj, the stochastic dispatch rule further reduces to: 


Rule # 

Pw 

Pi 

Pf 

applies if 

1 . 

1 

0 

0 

h > I 2 

0 < / < /2 

2 . 

I 2 

I-I 2 

0 


l2<l 

6 . 

1 

0 

0 


0<l<h 

7 . 

h 

0 

l-h 

h < I 2 


8 . 

h 

l-h 

h — h 


h<i 


The proof of Corollary 0] is provided in Appendix O 

For the sake of illustration, later on we will discuss examples for which the conditions of Corollary 
0] hold. Notice that, despite the fact that these conditions are quite restrictive, they prompt a 
simplified stochastic dispatch rule that still include cases for which the merit order is broken (see 
rules 7 and 8). 

To conclude this section, we would like to point out the case of symmetric balancing costs, that 
is, Cp — cf = Cp — cf = A. By Equation ([5|), we arrive at li = F“^(0.5), which is the median of 
the probability distribution that characterizes the stochastic power production. If, besides, this 
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probability distribution is symmetric, li is equal to its expected value. Under symmetric balancing 
costs and pp >Pw^ stochastic dispatch rule yields solutions that violate the merit order if 
C/ > 0.5cJ for the case in which cf> Cj. 

4. A Price-consistent Conventional Two-stage Market (ConvM-VB): The Role 
of Virtual Bidding 


The stochastic two-stage market StoM ensures maximum cost-efficiency and price-consistency. In 
contrast, it may violate the merit order and as a result, dispatch c onventional generating units in 


a loss-making position in the forward market ([Morales et al 


2014b|). In this section, we analyze an 


enhanced form of the conventional two-stage market ([I|)-(I1]) that, by construction, ensures price 
consistency and preserves the merit order. Mathematically, we can get price-consistent solutions 
out of ([I])-(I21) by freeing variable pw > 0, which represents the forward dispatch of the stochastic 
power capacity, while at the same time enforcing the additional constraint = Eq [A^(a;)]. As we 
show below, to do so, we need to pose the price-consistent conventional two-stage market as a 
complementarity problem. 

In practice, price-consistent solutions can be theoretically obtained from a conventional two- 
stage market where virtual bidding is allowed and exercised by a risk-neutral arbitrager that has 
perfect knowledge of the market price distribution (induced by the uncertain power supply). For 
this purpose, we solve the system of non-linear equations that results from concatenating the KKT 
conditions of the following optimization problems: 

Clearing of the forward market 


Minimize Cppp + Cjpi 

PF,PI,PW 

(12a) 

s.t. Pf + Pi + Pw + Pv = 1 ■ 

(12b) 

0 < Pw < Pw 

(12c) 

0<Pf<Pf 

(12d) 

0 < Pi < Pi 

(12e) 
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where pv is the optimal virtual bid coming from the arbitrager’s problem below. Notice that we 
assume that the bid price of the virtual bidder is zero. 

Clearing of the balancing market 


Minimize pf{u})cfp — Pp{uj)cp + vs{uj) (13a) 

pf{‘^),Pp0Apw0,^0 

s.t. pf{uj) — pf,{ui) + Apwioj) + Apviuj) + s{uj) = 0 : X^{u!) (13b) 

0<Pp{oo) <Pf-Pf (13c) 

0<pf{uj)<pF (13d) 

0 ^ Pw ^Pwi^) — 11 ^(^) ( 136 ) 

0<s{oj)<l (13f) 


where pw and pp are given by problem (fT^ and Apy{uj) is the amount of electricity resold, if 
positive, or repurchased, if negative, by the virtual bidder in the balancing market. 

Arbitrager’s problem 


Maximize pyX^ / Apy{uj)X^{uj)f{u})du} 
pvApvO Jq 

(14a) 

S.t. PyApy{oj) = 0 , 

(14b) 


where A-l is given as the Lagrange multiplier associated with constraint (112bl) and X^{ui) as the 
Lagrange multiplier associated with constraint (113bl) . 

We now characterize the solution to the system of KKT conditions associated with the convex 
optimization problems ffT^ . fflHll and (fTH) (the so-called short-run equilibrium solution). Notice 
that the forward dispatch associated with this solution exhibits a fundamental property, namely, 
it satisfies the merit order, while leading to price-consistency. 

Theorem 2 (A price-consistent conventional two-stage market). Consider the stylized 
power system described in Definition [H Consider also the equilibrium problem that results from 
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simultaneously enforcing the optimality conditions of problems m, d and m, where Cj <Cf- 
Define the constants 


Z2:=min I: [v - cf) F {I-pp) + cfF{l) > Cj (15) 

Z5:=min I: {v - cf) F {I-pp-pj) + cfF{l-pj) > Cp (16) 

Iq := min I : (y - cf) F {I -pp - pj) + cfF{l - pj) >Cf (17) 

and the function h ■ K>o —)• M>o, 


l7il) = F 


-1 


Cf- {v - cf) F {I - Pj - Pp) - CpF {I - Pj) 


Cp Cr 


Note that, by construction, h Fie- The equilibrium solution is given by: 


Rule # 

Pw+Pv 

Pi 

Pf 

applies if 

1 . 

1 

0 

0 

0<l<l2 

2 . 

h 

l-h 

0 

h < 1 ^l2 +Pj 

3. 

i-Pi 

Pi 

0 

I2 pj < / ^ /s 

4 - 

h{i) 

Pi 

l-h(l)-Pj 

I5 <l <Iq 

5. 

i-pp-Pi 

Pi 

Pf 

le < 1 


in all cases with pw < Pw ■ 


(18) 


A proof for this theorem can be found in Appendix IdI 

We can now use this theorem to provide conditions under which the conventional two-stage 
market with virtual bidding results in maximum cost-efficiency. 


Corollary 5 (Merit-order dispatch solution with virtual bidding). Consider the power 
system described in Definition [Jl where, in addition, we have that cj <Cf- The expected system 
operating cost associated with the price-consistent merit-order dispatch solution is equal to that of 
the stochastic dispatch solution in any of the following cases: 

1. Ii>l 2 and 0<l< min {pj + li,h); 

2. > A and <l <pj + li provided that pp > p^^ (capacity adequate power system). 

Cp—Cp Cp 

3. I ^ ih,Pw Fpj +Pp) if cf = Cf (free upward regulation). 
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4. 0<l<Pi+Pw or l>p^r+Pi+ Pf! 'ifcf = CF andpF>Pw (i-e., free upward regulation in a 
capacity adequate power system). 

5. CF = cf = cf (free upward and downward regulation). 

Furthermore, if Cf = Cp (free downward regulation), cf — Cf > 0, and F~^ (0) = 0, then the 
price-consistent conventional two-stage market does not deliver maximum cost-efficiency except for 
l>max{l^,li-\-pj) when l 3 >pp, and l>max{le,pp+pj) otherwise. 

The proof of Corollary [5] is provided in Appendix |El 

It is apparent that enforcing price consistency, through the introduction of virtual bidding, 
makes the dispatch solution of the conventional two-stage market substantially more intricate. Most 
importantly, virtual bidding enhances the efficiency of the conventional two-stage market to such an 
extent that the cost-efficiency gap with respect to the stochastic dispatch solution is closed under 
certain conditions. In essence, the ability of virtual bidding to close this gap is mostly determined 
by the relation between the characteristic system constants li and I 2 . On the one hand, the value 
of li (which is the solution to a news-vendor-type of problem as remarked in Appendix [B]) is driven 
by the relation between the cost of providing downward regulation versus the cost of providing 
upward regulation. On the other, the value of I 2 weighs the cost of dispatching the inflexible power 
capacity against the cost of dispatching the stochastic power capacity. The former action incurs a 
marginal cost of C/, while the latter entails a probable marginal cost of cf, or v in the worst-case 
scenario, given that any potential shortage of stochastic power production will have to be covered 
with upward regulation from the flexible generation technology and/or load curtailment. 

Thus, the effectiveness of virtual bidding to increase the cost-efficiency of the conventional two- 
stage market is greater for those power systems with a fleet of cheap inflexible power plants and 
comparatively costly means for downward regulation {li > I 2 ). In such a case, virtual bidding can 
even drive the conventional two-stage market to maximum cost-efficiency for a range of system 
loads. 

As previously mentioned, the conventional two-stage market with virtual bidding yields dis¬ 
patch solutions that respect the merit order, while resulting in price-consistency. This leads to an 
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interesting conclusion. Even under those conditions for which the merit order is preserved in the 
stochastic dispatch rule {li > 12 ), a price-consistent conventional two-stage market does not nec¬ 
essarily deliver maximum cost-efficiency. This happens, for example, in the load range h <l <Iq. 
This implies that a price-consistent merit-order market solution is not, in general, the merit-order 
market solution that minimizes the expected system operating cost. This will become more clear in 
the illustrative example below and when we introduce, later on, the conventional two-stage market 
with centralized dispatch of stochastic power production. 


4.1. Example 1 


In this example we consider the stylized power system described in Dehnition [T] with v = 
$1000/MWh. The uncertain power supply comes from a wind power farm whose capacity factor is 


assumed to follow a Beta distribution. The mean (k) and standard deviat ion (a) of this distri 


3ution 


are linked together through the empirical relationship ()19l) provided in 


Fabbrietal 


((20051). The 


shape parameters a and /3 of the Beta distribution modeling the wind power capacity factor are, 
consequently, computed according to ([^ . 


0-= 0.01837-^0.20355-K. 


(1 — k) ■ K - K 
a = --- K,, 




fi = a 


1 — K 


K 


(19) 

( 20 ) 


The expected wind power production is hence given by the product of the predicted wind power 
capacity factor k and the wind power capacity pyy. 

We now examine and compare the stochastic two-stage market and the conventional two-stage 
market with and without virtual bidding. For this purpose, we consider the three illustrative cases 
collated in Table [K and discussed below. Table [2] provides the dispatch solution {pi,Pf,Pw), or 
{Pi,Pf,Pw +Pv), and the expected total operating cost prompted by the efficient two-stage market 
and the conventional two-stage market with and without virtual bidding. This table also includes 
the forward price A-^ and the expected real-time price under each of the markets, and the 
incremental running cost in percentage with respect to the most cost-efficient dispatch solution. 
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Pi 

Pf 

Pw 

Cl 

Cp 

4 

Cp 

1 


Case a) 

500 

500 

100 

30 

35 

35 

30 

250 

0.5 

Case b) 

500 

500 

100 

30 

35 

40 

30 

250 

0.5 

Case c) 

100 

50 

100 

30 

35 

35 

30 

170 

0.5 


Table 1 Example 1: Data. System load and capacities are given in MW and marginal costs in $/MWh. 

which is the one provided by the stochastic two-stage market, logically. In the sequel, we will refer 
to this incremental cost simply as efficiency gap. 

• Case a) considers a power system with asymmetric balancing costs, in which upward regu¬ 
lation is less costly than downward regulation. More specifically, this instance corresponds to a 
capacity adequate power system with free upward regulation. Hence, the results in Corollary [2] 
apply and, accordingly, the stochastic two-stage market provides a solution that respects the merit 
order, because the capacity of the more expensive flexible power generation is reserved to deploy 
upward regulation in real time, when needed. Naturally, the conventional two-stage market always 
yields, by construction, a dispatch solution that satisfies the merit order. However, in the case that 
virtual bidding is not allowed, the conventional two-stage market clears the expected wind power 
production, which is, in general, a suboptimal dispatch decision in terms of market efficiency and 
as such, results in an efficiency gap greater than zero. On the contrary, the introduction of virtual 
bidding not only induces a price-consistent market solution, but also fully closes the efficiency gap. 
This result is consistent with Claim 4 in Corollary O since 0 <l <pj +Py^/, with I = 250 MW, 
Pj = 500 MW, and = 100 MW. 

• Case b) represents a capacity adequate power system identical to that of Case a), except for 
the fact that the provision of upward regulation is now costly, that is, c'^> Cp. In particular, we 
have that Cp — Cp = Cp — Cp = $5/MWh (symmetric balancing costs). Besides, since I < pj, the 
second result in Corollary [i] applies. Lastly, because C/ > 0.5cJ, the stochastic two-stage market 
provides a dispatch solution that breaks the merit order. Evidently, the conventional two-stage 
market cannot replicate such a dispatch, irrespective of whether virtual bidding is permitted or 
not. Nevertheless, virtual bidding reduces the efficiency gap, while ensuring price consistency. 
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Pw{+Pv) 

Pi 

Pf 


V 

cost 

gap 


StoM 

63.5 

186.5 

0 

30 

30 

6095 

- 

Case a) 

ConvM 

50 

200 

0 

30 

17.5 

6170 

1.2% 


ConvM-VB 

63.5 

186.5 

0 

30 

30 

6095 

0% 


StoM 

50.5 

188.5 

11 

30 

30 

6139 

- 

Case b) 

ConvM 

50 

200 

0 

30 

20 

6195 

0.9% 


ConvM-VB 

58.5 

191.5 

0 

30 

30 

6154 

0.2% 


StoM 

70 

100 

0 

37.09 

37.09 

3717 

- 

Case c) 

ConvM 

50 

100 

20 

35 

34.83 

3740 

0.6% 


ConvM-VB 

51.5 

100 

18.5 

35 

35 

3737 

0.5% 


Table 2 Example 1: Results. Power dispatch values are given in MW, prices in $/MWH, cost in $/h, and 


incremental cost (efficiency gap) in percentage. 

• Case c) is also similar to Case a), but with reduced generating capacities and system load. In 
fact, this instance corresponds to a capacity inadequate power system with free upward regulation. 
Consequently, results from Corollary[2]hold and the stochastic two-stage market prompts a dispatch 
solution that does comply with the merit order. However, as opposed to what happens in Case a), 
not even virtual bidding is able to close the efficiency gap under the conventional two-stage market 
in this instance. This shows that a price-consistent merit-order-based market does not necessarily 
deliver maximum cost-efficiency. 

To conclude this example, we investigate how the mean capacity factor k and the forecast horizon 
impacts the efficiency gap. To this aim, we consider Case d), for which data is provided in Table O 
Note that this case corresponds to a capacity adequate power system with asymmetric balancing 
costs. To be more precise, the provision of downward regulation in this system does not entail 
any extra operating cost since Cp = Cp = $35/MWh. Therefore, results from Corollary [H which 
in general render a dispatch solution that violates the merit order, apply. Accordingly, we should 
expect that the dispatch solutions induced by the conventional two-stage market, both with and 
without virtual bidding, feature a nonzero efficiency gap in this case. 
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Pi 

Pf 

Pw 

Cl 

Cp 


1 

Case d) 500 

500 

100 

30 

35 

40 35 

250 


Table 3 Example 1: Data for Case d). 


This is confirmed by the plots in Figure [H which illustrate such an efficiency gap in percentage 
(on the y-axis) for different values of k (on the x-axis) and forecast horizons (i.e., time spans in 
between the clearings of the forward and the real-time markets), namely, 1, 24 and 48 hours. Each 
forecast horizon is associated with a differe nt empirical r elatio nship between a and k in the form 


of (f^Oll . all of which have been taken from 


Fabbri et al, 


()2005ll . The continuous and dashed lines 


correspond to the cases with and without virtual bidding, respectively. 

Two observations are in order: First, the efficiency gap increases as the time distance in between 
the closures of the forward and the real-time market augments. This is hardly a surprising result 
related to the fact that uncertainty in wind power production grows as does the forecast lead 
time. Second, virtual bidding, or more generally, a price-consistent conventional two-stage market, 
can substantially reduce the efficiency gap as compared to the case of a conventional market that 
does not ensure price consistency, particularly in those situations where the contribution of the 
stochastic power production is expected to be important. 

Finally, it is worth noting that the efficiency gap decreases for high values of k in the case where 
virtual bidding is allowed. This is a direct consequence of the fact that, as the value of k grows, 
the resulting wind power distribution becomes more and more skewed towards high wind power 
production values. Intuitively speaking, the probability mass concentrates more and more towards 
values of wind power production closer to the wind farm capacity p-^ as k increases, a phenomenon 
that virtual bidding captures and takes advantage of. 

5. A Conventional Two-stage Market with Centralized Dispatch of the 
Stochastic Power Production (ConvM-CD) 


In the sequel, we analyze a variant of the conventional two-stage market dB)-®, whereby the 
forward dispatch pw of stochastic power capacity is centralized and determined by a non-profit and 
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Figure 1 Efficiency gap for different expected wind power capacity factors and forecast horizons 

all-knowing entity that seeks to minimize the expected total system operating cost. A transmission 
system operator, for example, could take on this role. The mathematical model t hat we present 


next to simulate this market organization is inspired from the one introduced in 


Morales et al. 


(j2014b[ l and takes the form of the following bilevel linear programming problem. 




Minimize 


CiPi + CfPf+ / {cpPp{u;)-CpPf{u;) + vs{uj)) f{uj)duj 


s.t. Pp{uj) — pp{uj) -h Apyy{uj) -|- s(w) = 0 : A^(w), Vw € O 
0 <Pp{uj) <pp — Pf, Vcj G n 
0<Pf{uj)<Pf, Vw € n 
0<Pw<Pw^ Vw G 0 
0 < pw + Apw{u}) < W(to), Vo; G n 
0<s(cu)<l, Vo; G n 


Minimize CfXf + CjXj 

Xp,Xj 


{Pf,Pi) Garg< 


s.t. Xf + Xj +pw = l ■ 
0<Xf<Pf 

0 < Xj <Pi 


. 


(21a) 

(21b) 

(21c) 

(21d) 

(21e) 

(21f) 

(21g) 

(21h) 

(21i) 

(21j) 

(21k) 
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Essentially, the lower-level problem (j21hp - (|21kp models the clearing of the forward market as 
a function of the amount pw of stochastic power production that is dispatched. This amount is 
determined by the upper-level problem (I21al) -(21g) with a view to minimizing the expected total 
system operating cost ()21al) . To this aim, the upper-level problem explicitly anticipates, through 
()21bp -( PT^ , the projected re-dispatch actions that will need to be undertaken for any possible 
realization oj of the stochastic power production. What is most important, though, is that the 
bilevel linear model (12 ip yields, by construction, a market solution that, while preserving the merit 
order (due to the enforcement of (121hp - (|21kp i. results in maximum cost-efficiency. 

The following theorem characterizes the solution to the bilevel linear programming problem (1211) . 


Theorem 3 (A conventional two-stage market with centralized dispatch of pw)- 

Consider the power system described in Definition{^with Ci <Cf- The solution to the bilevel linear 
programming problem m is given by: 


Rule # 

Pw 

Pi 

Pf 

applies if 

1. 

1 

0 

0 


0<l<l2 

2. 

h 

l-h 

0 


I 2 < 1 <Pj +12 

3. 

i-Pi 

Pi 

0 

h > h 

Pi + I 2 <l <Pi + li 

1 

h 

Pi 

l-Pj-h 


Pi + li <l <Pi +Pf +h 

5. 

I-Pi-Pf 

Pi 

Pf 


Pi +Pf +li<l 

6. 

1 

0 

0 


0 <l < min(/ 2 , (s) 

1. 

I 2 

I-I 2 

0 

li < I 2 

I 2 ^ 1 ^ Is 

8. 

h 

Pi 

l—pj — h 


Is ^ 1 pj + Pf “1“ 

9. 

1-Pi-Pf 

Pi 

Pf 


max{ls,Pi +Pp + I 1 ) < 1 


where constant /§ is defined when li < I 2 as: 


pmin{x-pj ,li+pp) 

Is '■= X : h +pj < X <l 2 +Pi and Cl {pj— m.a.x{x — l2,0)) — / {A{s) — B{s))ds— 

J minix—vr.h ) 


min{x—pj,li) 


rx-pi 


(A(s) — C{s))ds — 


^min(rc,Z2) 


A{s)ds = 0, 


( 22 ) 


' mm{x—pj fifi-pp) 


-Pl 


with 


(23) 
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A{s)\={v-c^)F{s-pp)-\-c%F{s) (24) 

B{s) := Cp + {v-cf) F{s-pp) - Cp (1 -F(s)) (25) 

C{s) := [v - Cp) F {s-pp)+ CpF (s) (26) 

A proof of this theorem is provided in Appendix [Fj 

As we did for the case of the conventional two-stage market with virtual bidding, we can now use 
Theorem [3] to identify conditions under which a conventional two-stage market with centralized 
dispatch of the stochastic power production delivers maximum cost-efficiency. 

Corollary 6 (Cost-efRcient merit-order dispatch solution). Consider the power system 
described in Definition{^ where, in addition, we have that Ci < Cp. The expected system operating 
cost associated with the merit-order dispatch solution provided by the bilevel linear programming 
problem m is equal to that of the stochastic dispatch solution in any of the following cases: 

1. Zi > Z 2 / 

2. Whenever the stochastic dispatch solution satisfies the merit order; 

3. Whenever the stochastic dispatch solution is such that pj =pj; 

4. cf = cj = cf (free upward and downward regulation). 

The proof of Corollary [6] is provided in Appendix [Gl 

Interestingly, the previous results show that all the stochastic dispatch solutions that do not 
violate the merit order could be recovered, in principle, from a conventional two-stage market, 
if we let an all-knowing and social-welfare-maximizer entity (for example, a TSO) control the 
dispatch of the stochastic power production. This implies that a TSO could close the efhciency 
gap between the stochastic and the conventional two-stage markets in many cases where virtual 
bidding fails to do so. On the other hand, neither a TSO nor virtual bidding are able to close such 
a gap in all those cases where maximum cost-efficiency requires breaking the merit order. In those 
cases, though, a TSO would be able to reduce the efficiency loss further than virtual bidding. In 
this vein, when compared with the results provided in Theorems [T] and O the dispatch solution 
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induced by a conventional two-stage market with virtual bidding reveals that price-consistency 
does not necessarily implies maximum cost-efficiency. But it is more interesting yet to notice that 
this conclusion also means that the conventional two-stage market with centralized dispatch of 
the stochastic power production can deliver dispatch solutions that are price inconsistent. Indeed, 
Theorem [2] tells us that this may happen in the range of loads in between and Iq. We illustrate 
one of these cases in the example of Section 15.11 

We conclude with the corollary below, which states that the price-consistent and the cost-efficient 
merit-order dispatch solutions are equivalent when the inflexible generating capacity is sufficiently 
large. 

Corollary 7 (Merit-order dispatch solution with large inflexible power capacity). 

Consider the stylized power system described in Definitions^ with Pj>l. Then, the solutions of the 
bilevel model m and the complementarity model consisting of the KKT conditions of problems 
(1121), dH and m are the same. The dispatch rule in that case is given by: 


Rule # 

Pw 

PI 

Pf 

applies if 

1. 

1 

0 

0 

0 < ^ < /2 

2. 

I 2 

I-I 2 

0 

h<i 


Proof. From Theorem O only rules 1 and 2 apply for pj > 1. From Theorem [3] with li > and 
Pj > I, only rules 1 and 2 apply. Furthermore, if li < I 2 and pj > I, then > I and consequently, 
only rules 6 and 7 in Theorem [3] apply. □ 

5.1. Example 2 

The main purpose of this example is to illustrate the key differences between the merit-order 
dispatch solutions provided by the two enhanced variants of the conventional two-stage market 
that we have examined in this paper, namely, that in which virtual bidding is allowed (ConvM-VB) 
and that in which the stochastic power production is centrally dispatched (ConvM-CD). To this 
aim, we consider Case e), for which data is provided in Table 01 Note that this new case is identical 
to Case d), but for a different system load of 155 MW. Results for the several markets analyzed are 
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collated in Table [5l These results consist of dispatched quantities (pw^Pi^Pf) or {pw +PviPfPf), 
as appropriate; forward price (A-f); expected real-time price (A**); expected total system operating 
cost (cost), and efficiency gap. 

Case e) corresponds to a capacity inadequate power system with asymmetric balancing costs. In 
particular, the provision of downward regulation in this system is free in the sense that Cp = Cp. 
Therefore, the stochastic dispatch rule is governed by the results of Corollary [H which generally 
prompt dispatch solutions that break the merit order, as is the case here (see row labeled “StoM” in 
Table [5|). Hence, none of the enhanced variants of the conventional two-stage market we consider, 
i.e., ConvM-VB and ConvM-CD, are able to close the efficiency gap (see column tagged as “gap” 
in Table [5]). 

The stochastic two-stage market delivers minimum expected system operating cost at the same 
time that it guarantees price consistency. However, in order to produce a price-consistent and max¬ 
imum cost-efficient dispatch solution, it has to violate the merit order. This has impo rtant negative 


implic ations towards the actual implementation of this market-clearing procedure ([Morales et al. 


2014b|). As an example, it is easy to infer from the results in Table [5] that the flexible power 


capacity is dispatched in a loss-making position under the stochastic two-stage market, because 
the forward price equals $30/MWh, whereas the marginal production cost of the flexible power 
generation technology is $35/MWh. 

On the other hand, the conventional two-stage market with virtual bidding (ConvM-VB) pro¬ 
duces a price-consistent dispatch solution that respects the merit order, but that is not the best 
in terms of cost efficiency. Indeed, the market solution provided by ConvM-VB is even less cost- 
efficient than the one delivered by the plain conventional two-stage market (ConvM). Although 
this should be regarded as a rare case, it shows that ensuring price consistency does not necessarily 
lead to higher market efhciency (understood as the minimization of system operating costs). 

Lastly, the conventional market with centralized dispatch of the stochastic power production 
(ConvM-CD) yields the most cost-efficient dispatch solution among those that comply with the 
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Pi 

Pf Pw 

Cl 

Cf cf 

Cf 

1 

K 

in $/MWh. 

Case e) 100 

50 100 

30 

35 40 

35 

155 

0.5 

Table 4 Example 2: Data. System load and capacities 

are given in MW and marginal costs 


Pw i+Pv) 

Pi 

Pf 

A/ 


cost 

gap 

StoM 

18.5 

92.5 

44 

30 

30 

3245 

- 

ConvM 

50 

100 

5 

35 

25.42 

3296 

1.6% 

Case e) 








ConvM-VB 

58.5 

96.5 

0 

30 

30 

3304 

1.8% 

ConvM-CD 

18.5 

100 

36.5 

35 

22.93 

3272 

0.8% 


Table 5 Example 2: Results. Power dispatch values are given in MW, prices in $/MWh, cost in $/h, and 


incremental cost (efficiency gap) in percentage. 

merit order. Note that ConvM-CD gets to reduce the efficiency gap by half. To do so, however, it 
must relinquish price-consistency. 

In short, preserving the merit order in forward electricity markets with uncertain supply implies 
giving up on cost-efficiency or price-consistency (or both). 

6. Conclusions 

The overall message that emerges from our analysis is that the concepts of merit order, cost- 
efficiency, and price-consistency are conflicting requirements that cannot be met together in an 
electricity market. Indeed, preserving the merit order generally comes at the expense of market 
efficiency, either in the form of price-inconsistent or cost-inefficient market solutions. 

We have reached this conclusion after examining four types of two-stage markets with uncertain 
supply, in which only one or two of the above requirements can be guaranteed in general. We 
have provided mathematical formulations for these four types of markets and explain how these 
mathematical abstractions translate into practice. 

Our study has also revealed general conditions under which preserving the merit order is most 
likely to jeopardize cost-efficiency and/or price-consistency. Interestingly, this seems to be the case 
of power systems that are capacity inadequate (broadly understood as a power system where the 
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complete lost of the uncertain supply cannot be covered with economical flexible power generation) 
and/or where the provision of upward regulation is costly, while the provision of downward regu¬ 
lation comparatively imposes little or no extra cost at all to the system. Likewise, our study has 
identified general conditions under which the merit order, price-consistency and cost-efficiency can 
indeed be met together. This is, for example, the case of power systems where the implementation 
of real-time adjustments do not entail opportunity costs. 

In an attempt to keep our investigation essentially analytical, we have built our market models 
on a stylized power system with inhnite transmission capacity. Consequently, a natural avenue for 
future research is to elucidate whether more realistic assumptions on the underlying power system 
can diminish the loss of efficiency caused by the merit order or even limit the cases of inefficiency 
to a few “degenerate” ones. To this end, we are most likely to abandon our analytical approach 
and make use of computational simulation instead. 

On a different front, our analysis has considered energy-only electricity markets. Therefore, 
another logical direction for future work is to investigate whether the consideration of other types 
of market mechanisms could reduce or even nullify the loss of efficiency associated with the preser¬ 
vation of the merit order. In this regard, we conjecture, based on our results, that appropriate 
markets for downward operating reserve might do the trick. 

Appendix A: Proof of Proposition [2] 

Proof. For given values of pp and pw, the optimal re-dispatch actions 
pf{uj),p'f{uj),Ap\Y(u;),s{Lo) are determined by solving the following optimization problem. 


Minimize / (vs{u) + cfpf.{uj) — CpPp{u})) f{uj)duj (27a) 

s.t. s{u})+pf{oj)—pf{oj) + Apyp(oj) = 0, Vo; G n (27b) 

0<Pf{u})<PF, Vw G n (27c) 

0 <pf{oj) <pp — pp, Vo; G 0 (27d) 

0 ^ Pw T Apw (w) ^ IF ((u), G n (27e) 

0 < s{u}) <Pf + Pwj Vo; G II (27f) 
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Depending on whether the realized stochastic generation W{oj) is higher or lower than the 
dispatched quantity pw, the re-dispatch rule for each scenario uj is given by: 


UW{lo) <Pw 


pf{u})= min {Pf-Pf,Pw -W{u})) 

< p-[u)) = 0 

^s{uj) = max(0,pn/ -W{u}) -pp +Pf) 


(28) 


liW{Lo)>Pw< 

Therefore, the second-stage expected cost can be computed as: 


P+{lo) = 0 


= min {pf,W{uj) -pw) 
s(a;) = 0 


(29) 


{ys{u) -h CpPp{u) - CpPp{u)) f{u)du = 

POO 

{vs{u;) + CpPp{uj)-CpPp{u;)) f{u;)duj+ / {vs{uj) + CpPp{u;) - CpPp{u;)) f{u;)du} = 


rPw 


'Pw 


rPw 


rPw 


vs{u})f{ui)doj+ / cj,pp{uj)f{uj)duj- / CpPp{u;)f{uj)doj = 

Jo J p\Y 


fPF+PW-PF 


= V 


{pw -W{u})-Pp+ Pf) fi‘^)duj + c 


fPF+PW-PF 


(Pf-Pf) f{^^)du}+ 


+ c] 


rPw 

■ 


{pw - W(uj)) f{u})du} - c} 


rPF+Pw 


{W{uj) - Pw) f {oj)duj - Cj 


PFfiuj)duj = 


PF+PW-PF 


= v{pF+Pw-PF)FiPF+Pw-PF)-'V 


•'PW 
rPF+Pw-PF 


’PF+PW 


W{ijj)f{uj)duj+ 


rPw 

+ cj. {pf - Pf - Pw) F {pf + Pw - Pf) + c%PwF (pw) -c^ 

J or+n 


+ Cf [Pf + Pw) F {pf +Pw) - CfPwF {pw) - CfPf - Cf 

rPF+PW-PF 


W {u})f{uj)duj+ 

PF+PW-PF 

pf+pw 

W {u})f{Lo)dLO = 


= V 


F{oj)doj + c 


rPw 


•'Pw 

PF+PW 


F{uj)doj + Cf 


F{uj)duj — CfPf 


(30) 


' PF+PW-PF 


'Pw 


where, in the last equality, we have used the integration-by-parts theorem, according to which 
Ixf ^fi^)ds = X 2 F (X 2 ) - XiF (xi) - F (s) ds. □ 


Appendix B: Proof of Theorem [T] 

The proof of Theorem [1] relies on the proposition presented and proven below. 


Proposition 4. Consider a power system as deseribed in Definition{l\in whieh the dispateh of the 
inflexible unit pj is given. Let p then denote the net load to be satisfied by the flexible and stochastic 
power generating units. The dispatch rule for the flexible and stochastic generation is given by: 
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Rule + 

Pw 

Pf 

applies if 

1. 

P 

0 

0 <p < h 

2. 

h 

p-h 

h<p<h +Pf 

3. 

P-Pf 

Pf 

h +Pf 


Likewise, the marginal production cost of the flexible-stochastic generation portfolio, denoted by 
c, writes as: 

{v-cf)F{p-pp) + cfF{p) if 0<p<li 

c= < CF-\-{v-cf)F{p-pp)-cf{l-F{p)) if li<p<li-\-pp (31) 

^{v-c~p)F{p-pp)-^cfF{p} if li-^pp<p 

Proof. Using Proposition^ the optimal forward dispatch of a portfolio of flexible and stochastic 
power generation is given as the solution to the following optimization problem: 


rP-PF 


Minimize z = Cppp + v 
s.t. 0<Pf<Pf-{iF) 


F{u)duj + cl 


fp-PF 


F{ijj)duj + Cj. 


'P-PF 


'P-PF 


F{uj)duj — CfPf (32a) 
(32b) 


where 7,7 are dual variables and pw =p — PF- Problem (|32l) is a convex optimization problem that 
satishes a Slater condition and therefore, the KKT conditions below are necessary and sufficient 


for optimality. 


dz _ 

- -7 + 7 = 0 

opF - 

0 < Pf -L 7 > 0 

0<{Pf-Pf) 


(33a) 

(33b) 

(33c) 


where 

Fi y 

— = {cf - Cf) - (c+ - cf) F (p-Pf) (34) 

OpF 

We now determine the optimal solution to problem fl32j) by exhaustively enumerating the points 
that satisfy the optimality conditions (|33p as follows: 
a) Pf = 0 


(|33cD —5-7 = 0 

F{p)<%^ 

OpF cf — Cp 


(i33bl)^o' = n 


(l33cD —>7 = 0 

(I33ap -5- = 0 F{p-pp) = ^ 

OpF Cf—Cp 


b) 0 <Pf <Pf 
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c) Pf=Pf 




i^^F{p-pp)>^^^ 


Cp Cl 


Since pw =p — PF, the solutions of cases a)-c) can be summarized as the following dispatch rule: 



p if 

Pw = Ih if 

P-Pp if 

0<p<li 

li K p K li F Pp 

h Fpp<p 


(35) 

where h = F ^ ). The marginal production cost of the flexible-stochastic generation port- 

\ Cp—Cp ) 

folio c is equal to ||, that is. 

c= < 

(u - cj,) F {p- pp) + c+F (^ 

Cf + {v - cf) F {p -pp) - Cp (1 
^ {v - Cp) F {p- Pp) -h CpF (p) 

if 

-Fip)) if 

if 

0 <p < Zi 
li<p<li Fpp 
liFpp <p 

(36) 


Note that the functions pw (p) and c (p) are increasing and continuous on p. 


□ 


Remark 1 (News-vendor solution). The characteristic con stant L in (15^ and (1551) can be 
interpreted as the solution to the classical news-vendor problem (jRaiffa and Schlaifeiil2000ll : 


where Qopt is the optimal stocking quantity of the news-vendor, F {■) is the cumulative distribution 
function of the demand to be satisfied, c is the variable production cost, and p and h correspond 
to the penalty cost of unsatisfied orders and the inventory holding cost, respectively. The analogy 
works, thus, as follows: The variable cost of the stochastic power production is zero, i.e., c = 0; p = 
Cf — Cp represents the penalty cost of dispatching less stochastic power capacity than its eventual 
power production, since the consequent power surplus is to be compensated for by a decrease in 
the flexible power generation (downward regulation); finally, h = Cp — cf provides the marginal 
cost of dispatching more stochastic power capacity than its eventual real-time power production, 
because the consequent generation deficit is to be covered with an increase in the flexible power 
generation (upward regulation). Therefore, 


Qopt ^ 


p — c 
p + h 


= F- 


CF-Cp 


Cf Cp F Cp Cp 


= F- 


Cp Cl 


^ =h. 


(38) 



















36 


Morales and Pineda: On the Inefficiency of the Merit Order in Forward Markets with Uncertain Supply 
Article submitted to Operations Research^ manuscript no. (Please, provide the manuscript number!) 


Remark 2 (Discrete probability distribution). Proposition H assumes a continuous prob¬ 
ability distribution for the uncertain electricity supply, which implies that objective function (I32ap 
is differentiable. If the uncertain supply is, in contrast, characterized by a discrete probability dis¬ 
tribution, the cumulative distribution function F {■) is stepwise and thus, objective function (I32ap 
becomes nondifferentiable. In that case, the sub-derivative of the objective function should be used 
instead to formulate the KKT optimality conditions as follows: 


dz- _ dz+ 

+ 7 — 7 < 0 < —-1- 7 — 7 


dpF 
0 <Pf -L 7 > 0 


dpF 

0<{Pf-Pf) 


where 


dJ+ 

dpF 

dz~ 

dpF 


= CF-Cp- (c^ - Cp) F{p- Pp) 

= cp-cp- {c%-cp) {F{p-pp)-f{p-pp)) 


(39a) 

(39b) 

(39c) 

(40a) 

(40b) 


Note that /(•) should be interpreted here as the probability mass function. We analyze next the 
different points satisfying the optimality conditions 
a) Pp = 0 


(M 

(I39ap 


'7 = 0 

dz+ 


— 7 > 0 


dz+ 

dpi 


>0 


F\ 


< 


Cf - Cf 


— 


Cp Cl 


b) 0<Pf<Pf 


(Ml 

(Ml - 


dpF - 

-A 7 = 0 

-7> 7 = 0 


■cF{p-pp)<^ — ^ + f{p-PF) 


Cp Cl 


Cp Cr 


C) Pf=Pf 


(I39ap 


'7 = 0 
dz 


dpi 


+ 7<0 ^ CF-Cp- (c+-Cp,) {F{p-pf)- f {p-Pf))+1<^ 

Cp Cp 


F{p-Pp) > 


Ctp Cz 


Since F ^ (•) is the generalized inverse distribution function, and li is defined as F ^ ), 

\ Cp—Cp, J 

the dispatch rule for a discrete probability distribution of the uncertain supply coincides with ([5^ . 
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Using Proposition m the proof of Theorem [T] proceeds as follows: 

Proof. This proof deals with the optimal dispatch of an inflexible power unit with capacity 
Pj > 0 and marginal cost Cj and the flexible-stochastic generation portfolio of Proposition [H For 
a total system load denoted by I, this optimal dispatch is determined by solving the following 
optimization problem: 

Minimize z = 

p/>0,p>0 

Pi+P = 1:t (41b) 

Pi<Pj\(l) (41c) 

Since the integral of an increasing function is a convex function, problem m is a convex 
optimization problem and therefore, the KKT conditions below are necessary and sufficient for 
optimality. 


r~ 

CiPi + / c{x)dx (41a) 

Jo 


0 < {cj — T + (j)) j- Pi > 0 

(42a) 

0 < (c{p) — t)±p>0 

(42b) 

0< (pj-p/) T(?i>0 

(42c) 

pi+p = l 

(42d) 


The solution to problem fl41j) is then obtained by exhaustively examining the points that satisfy 
the optimality conditions (l4^ as follows: 
a) Pj = 0 and p = 0 

(I42dp — )■ only feasible if ^ = 0 
(liSa ^d> = 0 
(I42ap -^T <cj 
(l4^ ^ r < c(0) 

h) Pi = 0 and p> 0 


()42dl) -^p = l 

(j42iP^(/) = 0 

(|42ap —>■ r = c(p) 


()42bp T <ci 


c{p) < Cj 
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c) Pi = Pi and p = 0 





Cl <T < c(0) 


d) Pi = Pi and p> 0 


(|42d[) p = l — Pi 


m 


• 0 > 0 

> Cl 


e) 0 < Pi < Pi and p = 0 


(I42bp — >• r = c{p) ^ 

(142d[l —)■ Pi = I 
(Ii2a^d) = 0 

fl42ap T = Cl 
(lM ^r<Hn'l 


Cl < c( 


Cl < c(0) 


f) 0 <pi <Pi and p>0 


Ci=C{ 


(I42dp -^pi +p = l 

(I42cp ^0 = 0 
(l42aD -^T = Cl ) 

(I42bp -^T = c{p) j 

Let p denote the value of p such that c(p) = C/. Note that the function c{p) is increasing and 
continuous on p and therefore, p exists provided that c(0) < Ci < c(oo). If C/ < c(0), then we set 
p = 0. Similarly, if c/ > c(oo), we assign the value oo to p. This way the solutions obtained in cases 
a)-f) above can be summarized in the following dispatch rule: 

I if 0 <l <p 

p=lp if p<l<p + Pi (43) 

I- Pi if p + Pi<l 

where the dispatch of the inflexible power capacity is pi = I — p. Note that the function c{p) is 
piecewise and therefore, p can take on the following values 

"O if c/<c(0) 

I 2 if c{0) < Cl <c{li) 

p=<h if c{li) <ci <c{li+pp) (44) 

li if c(Zi+p^) < C/< c(oo) 

^00 if c(oo) < Cl 
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where 

l2-.= l.{v- Cp) F{1- pp) + CpF {1) = Cl 
h\=l-.[v- cj) F{1- Pp) -Cp{l-F (1)) = ci-cp 
h:=l\{y- Cp) F{l-pp) + CpF (1) = cj 
Merging ([55)1 . ()551) and (|551) we obtain the following cases: 


Rule # 

Pw 

Pi 

Pf 

applies if 

1 . 

0 

1 

0 



0 ^l ^ Pi 

2 . 

l-Pi 

Pi 

0 

Cl < c( 0 ) 


Pi <l <h+Pi 

3. 

h 

Pi 

l-h-Pj 



h +Pi <l<h +Pi +Pf 

4. 

l-Pi-P f 

Pi 

Pf 



h +Pi +Pf < 1 

5. 

1 

0 

0 



0<l<l2 

6 . 

I 2 

l-h 

0 



I 2 ^ ^ ^ ^2 Pj 

7. 

l-Pi 

Pi 

0 

c( 0 ) <ci< c{h 

) 

I 2 + Pi ^ ^ -|- Pi 

8 . 

h 

Pi 

l-h-Pj 



h +Pi <l<h +Pi +Pf 

9. 

1-Pi-Pf 

Pi 

Pf 



h +Pi +Pf < 1 

10 . 

1 

0 

0 



0<l<h 

11 . 

h 

0 

l-h 



< ^3 

12 . 

h 

l-h 

h — h 

c{h) < Cl <c{l 

i+Pf) 

Z 3 < Z < I 3 +Pj 

13. 

h 

Pi 

l-h-P i 



h +Pi <l<h +Pf+Pi 

14. 

1-Pi-Pf 

Pi 

Pf 



h +Pf+Pi < ( 

15. 

1 

0 

0 



0<l<h 

16. 

h 

0 

l-h 



h <l <h +Pf 

17. 

1-Pf 

0 

Pf 

c{li +Pf) < c/ 

< c(oo) 

Pf ^ ^ ^ ^4 

18. 

I4 Pp’ 

l-h 

Pf 



I4 1 ^ I4 Pi 

19. 

l-Pi-P f 

Pi 

Pf 



I4 + pj ^ 1 

20 . 

1 

0 

0 



0<l<h 

21 . 

h 

0 

l-h 

c(oo) < Cl 


h<l<h +Pf 

22 . 

1-Pf 

0 

Pf 



h ~\~Pf — ^ 


Next the different dispatch rules above are recast as a function of the system characteristics 
constants only. To do so, we note that if hlzli do not exit, their values are assigned to 

infinity. 
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Rules 1-4 only apply if the marginal cost of the stochastic-flexible portfolio for p = 0 is higher 
than the marginal cost of the inflexible power unit C/. Note also that the condition cj is 

equivalent to li > I 2 and thus we can jointly reformulate rules 1-9 as follows: 


Pw 

Pi 

Pf 

applies if 

1 

0 

0 


0<l<l2 

h 

1 —12 0 


I 2 ^ 1 ^ Pj I 2 

l-Pi 

Pi 

0 

l\ > I 2 

Pi + l2<l<Pj + h 

h 

Pi 

l-h-Pj 


Pi + h < 1 <Pp +Pj + h 

1-Pf-Pi 

Pi 

Pf 


^ ^ Pf ~^Pi + b 


Rules 10,15 and 20 can be easily merged as 


Pw Pi Pf 

applies if 

1 0 0 

b < h 

0 <l<li 


The condition c{li) <ci< c{li +Pp) can be equivalently formulated as li < I 2 and Z 3 < li +Pp., 
which allows expressing rules 11-14 as: 


Pw Pi Pf 

applies if 

li 0 1 — li 

li 1 — I3 h — li 

h Pi l-h-pi 

1 -Pi-Pf Pi Pf 

h < h 

I3 <Ii+Pf 

< /s 

^3 < / ^ Pj + I3 

h+Pi <l<Pi-\-PF + h 

l>PJ-\-pp-\-ll 


Rules 16 and 21 can also be merged as: 


Pw Pi Pf 

applies if 

0 1 — /i 

h < h 

^3 > ^1 + Pf 

/i < Z < pp + / i 


Likewise, rules 17-19 and 22 can be rewritten as: 


Pw Pi Pf 

applies if 

1-Pf 0 Pf 

b Pp ^ b Pp 

1-Pf-Pi Pi Pf 

l-Pp 0 Pp 

b < b 

I 3 >li+Pp 

3b 

(ci < v) 

Pf “1“ Zi ^ Z ^ Z 4 

Z 4 ^ Z ^ Z 4 “h Pj 

Z > Z 4 + pj 

$l4 

l>Pp + h 


The final dispatch rule can be thus summarized as: 















Morales and Pineda: On the Inefficiency of the Merit Order in Forward Markets with Uncertain Supply 
Article submitted to Operations Research-, manuscript no. (Please, provide the manuscript number!) 


41 


Rule # 

Pw 

Pi 

Pp 

applies if 

1. 

1 

0 

0 


0<l<l2 



2. 

I 2 

l-h 

0 


^2 < ^ ^ P/ H" ^2 



3. 

l-Pi 

Pi 

0 

h > h 

Pi + h <l <Pi + h 



4. 

h 

Pi 

l — h—pj 


Pi + h <l <Pp +Pi +h 


5. 

l-Pp-Pi 

Pi 

Pp 


1 > Pp + Pi + h 



6. 

1 

0 

0 


0<l<h 

7. 

h 

0 

l-h 



VI 

VI 


8. 

h 

l-h 

h — h 



h <■ 1 ^Pi + h 






h <h+Pp 



9. 

h 

Pi 

l-h-P i 



h +Pi < 1 <Pi +Pp + h 

10. 

1 — Pi ~Pf 

Pi 

Pp 

h < h 


1 >Pi +Pp + h 

11. 

h 

0 

l-h 



h <l <Pp + h 

12. 

1 — Pp 

0 

Pp 



3U 

Pp + h<l<U 

13. 

U -Pp 

l-U 

Pp 


h > h +Pp 

{ci < v) 

^4 ^ ^4 Pj 

14. 

1 

1 

Pi 

Pp 




1 'f> 1 4 pj 

15 . 

i-Pp 

0 

Pp 



$U 

1 ^ Pf -\-li 


□ 

Remark 3 (Discrete probability distribution). The same sub-differential analysis used in 
Proposition m can be applied here so that the dispatch rule above is also valid for the case in which 
the uncertain power supply is modeled by a discrete probability distribution. In such a case, it 
suffices to redefine the constants as: 

I2 := min l\{y- cf) F {I-pp) + c^F{l) > C/ 

^3 := min I ■. Cp + {y - cf) F {I -pp) - Cp {I - F{ 1 )) > Cj 
^4 := min. I: (v — Cp) F {I—pp) + CpF{l) > Cj 

i>ii+pp ' ^ 

Appendix C: Proof of Corollary [4] 

Proof. First, we show that, under the conditions stated in this corollary, < li +pp, and there¬ 
fore, rules 11-15 of the stochastic dispatch solution do not apply. For this purpose, consider expres¬ 
sion m, which defines constant I 3 , and note that Cp + (^v — cf) F {pp — pp) — Cp{l — F{pp)) = 
Cp + {v — cf)F (0), because F (pp) = 1 given that pp > p^/. Furthermore, it holds that -|- (u — 
cf)F (0) > C/, since Cp > C/ and {v — cf)F (0) > 0. Consequently, we have that h<Pp < h +Pp. 
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Now we prove that ^3 = max {h^A) with A = F~^ (l — j is optimal for the minimization 
problem (ITOl) . By construction, < p^y. Note that Cp + (^v — cf) F (A — pp) — Cp{l — F{A)) — Cj = 
[v — cfj F {A — pp) > 0 and that F (x) = 0 for all x < 0. Hence, either l^ = A \i A>li or Z3 = li 
otherwise. 

We can proceed analogously to show that I 2 = F~^ (^) ’ denote constant F~^ (?) 
construction 0 < H <Pyy Consider minimization problem Q. It holds that (u — cJ) F {B —pp) + 
CpF (B) — cj = (v — cJ) F {B —pp) > 0. Since B — pp < 0 and F (x) = 0 for all x < 0, I 2 = B is 
optimal for dH). 

Finally, the result corresponding to I <pj trivially follows by isolating those ranges of system 
load for which I must be lower than pj. □ 

Appendix D: Proof of Theorem [2] 

Proof. Let us start identifying the optimal solution to the forward dispatch model (1121) for a 
given value of the virtual bid volume pv The merit-order dispatch of the inflexible, flexible and 
stochastic power capacity is provided in the table below. Note that the case I >pw + Pv +Vi +Vf 
makes problem (|12p infeasible. 


case 

PF 

Pi 

Pw 


applies if 

a) 

0 

0 

l-pv 

0 

1 <Pw +Pv 

b) 

0 

0 

Pw 

[0,cd 

l=Pw +Pv 

c) 

0 

l-pw -pv 

Pw 

Cl 

Pw +Pv <l<pw +PV +Pi 

d) 

0 

Pi 

Pw 

[ci,Cf] 

l=Pw +Pv +Pi 

e) 

l-pw-pv-Pi 

Pi 

Pw 

Cf 

Pw +Pv +Pi <l <Pw+Pv+Pi+Pf 

f) 

Pf 

Pi 

Pw 

[cf,oo] 

1=Pw+Pv+Pi+Pf 


Observe that in cases a), b) or c), the flexible generation technology is not dispatched and 
therefore, the balancing market-clearing problem (I13p simplihes to: 


Minimize cfp'f{u) + vs{u}) 

(45a) 

s.t. p'fiui) + /S.pw{pj) + + Apy{u}) = 0 : A^(a;) 

(45b) 

0<pf{u}) <pp 

(45c) 

0 < Pw + ^Pw (<^) < kF (cj) 

(45d) 

0 < s(cu) < 1 

(45e) 


whose solution boils down to a merit-order-based dispatch, that is. 
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Pf(^) 

Pf{^) 

Apvr(a;) 

s{uj) 

\\u^) 

applies if 

0 

0 

Pv 

0 

0 

Pw+Pv < W(uj) 

0 

0 

Pv 

0 

[0.C+] 

W{oj)=pw+Pv 

Pw +PV — W{uj) 

0 

W{uj)-pw 

0 

4 

Pw+Pv -Pf <W(w) <pw -I-Pv 

pw +pv — W{uj) 

0 

W (oj) — pw 

0 

[4.'y] 

W{oj)=pw +pv —Pf 

Pf 

0 

W{uj)-pw 

Pw +Pv -W{uj)-Pf 

V 

W(u>) <Pw +Pv -Pf 


The expected balancing price for these cases is computed as: 

r rPw+Pv-PF rPw+Pv 

/ \^{uj)f{ijj)doj= / vf{u})du}+ / Cpf{uj)duj = vF{pw+Pv—PF) + 

d do dpyy+pY—pp 

+ Cf {F {pw + Pv) - F {pw + Pv -Pf)) = {v - cf) F {pw + Pv-P f) + CfF{pw+Pv) 

In cases d), e) and f) of the forward dispatch, pw =Pw and therefore, the balancing market¬ 
clearing problem (II3D becomes: 


Minimize cjpj(w) — CpPp{(jj) + vs{uj) (46a) 

p%{‘^),Pp{‘^)Apw{‘^)^^(.‘^) 

s.t. p'^{uj)—p~p{uj) +/S.pY/{oj) + s{uj) +/S.pv{oj) = ^'- (46b) 

^<Pf{oj)<Pf (46c) 

Q<pt{oj) <Pf-Pf (46d) 

0 <pw + Apw(w) < 4k(w) (46e) 

0 < s{oj) < I (46f) 

whose solution is given by: 


whose solution is given by: 


Pf4) 

Pf{^) 

Apw (a;) 

s(uj) 

X\u;) 

applies if 

0 

1 

1 

1 

1 

1 

0 

0 

l-Pj<W{uj) 

0 

1 

1 

1 

1 

1 

0 

[0, Cp] 

W{uj) = l-pj 

0 

W{ijj) -pw -pv 

W (w) —pw 

0 

Cp 

Pw + Pv <W{uj) <l - Pi 

0 

0 

Pv 

0 

[Cf 1 +f] 

W(uj) =pw +Pv 

Pw +Pv - W(uj) 

0 

W (uj) —pw 

0 

Cp 

1 -Pi - Pf <W{lu) <pw +Pv 

Pw +Pv - l+Pi +Pf 

0 

1 -Pi -Pf -Pw 

0 

[ct^] 

W{uj) = 1-Pi-Pf 

Pf - 1 +Pw +Pv +Pi 

0 

W{oj) —pw 

1-Pf-Pi-W{uj) 

V 

W{lo)<1-Pj-Pf 
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Hence, the expected balancing price in cases d), e) and f) is calculated as: 

r pl-Pi-PF rPw+Pv fl-Pi 

/ X^{u})f{uj)duj= / vf{uj)du}+ / Cpf{oj)duj+ / Cpf{oj)duj+ / Of{u})du} = 

J Q Jo dl—pj—pp Jpw+Pv dl—pj 

= vF{l-pj- pp) + c% {F {pw +Pv)-F{l-pj- pp)) + Cp {F {I -pj)-F {pw +Pv)) = 

= {v - cj,) F {I -Pj -pp) + {cj, - Cp) F {pw + Pv) + CpF {I -Pj) 

Finally, the KKT conditions of the arbitrager’s problem (fn|) imply that: 
d 


dp 


V 


(^X^Pv-pv j 


)^{uj)f{ijj)doj ) =0 X^= X!’{u)f{uj)duj 


(47) 


That is, the strategy of the arbitrager is to place a zero-price virtual bid pv in the forward 
market and to repurchase or resell the same amount in the balancing market so that the forward 
price equals the expected balancing price. 

The solution to the short-run equilibrium problem must be a point that is simultaneously optimal 
for the arbitrager’s problem and the clearing problems of the forward and the balancing markets. 
Therefore, in the short-run equilibrium solution the forward price A-f must be equal to the expected 
balancing price J^X^{io)f{uj)du}, which we denote by X’’. Below we analyze the solution to the 
short-run equilibrium problem for each of the cases a) to f) included in the first table of this proof. 

a) 

A^ = 0 'I 

W = 0 


b) 


c) 


X^ = {v- cj,)F {I - Pp) + CpF (1) ^ 
0<Xf < Cl 

F^={v-ci)F{l-pp) + c+F{l)^ 


0 <1<L 


■ Pw +Pv — h 


X^ = Cj 

= {v - c-^) F {pw + Pv -Pp) + CpF {pw +Pv). 

Pw+Pv <l <Pw+Pv+Pi. 


I 2 F I I 2 F Pj 


d) 


Cl < X^ < Cp 

X^ = {v - cj,) F {I -pj -pp) + cj,F {I -pj) ^ 


I2 F I ~ Pi F h 


e) 


xf = 


Cf 


Xb = (^v - cf) F {I - Pj -pp) [c'^-Cp)F{pw+Pv) + CpF{l-pj)^ 

cp-{v- cj) F{l-pj-pp)- CpF {I - pj) 


F'{Pw+Pv) = 


Cf Cr 


Pw+Pv+Pi <1<Pw+Pv+Pi+Pf 
I -Pi -Pf <Pw + Pv <l -Pi 
F{l-Pi-P f)<F{pw+Pv) <F{l-pj) 


> < I — Pi < Iq 
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f) 


where 


Cp < A-^ 

A'' = {v - cf) F {I -pj -pp) + CpF {I -pj) 


I 2 ■■= min 1: [v - Cp) F {I -pp) + CpF{l) > Cj 

I5 := min I : [v - Cp) F (l-pp-pj) +CpF{l-pj) > cp 

le ■■= min I: [v - Cp) F {I -pp - pj) + CpF{l - pj) > Cp 

Note that and Iq are defined using minimization problems to account for both continuous and 
discrete probability distributions for the uncertain power supply. Therefore, the dispatch rule can 
be summarized as: 


Rule # 

Pw+Pv 

Pi 

Pf 

applies if 

1. 

1 

0 

0 

0<l<l2 

2. 

h 

I-I 2 

0 

I 2 <l <l2 +Pj 

3. 

i-Pi 

Pi 

0 

I 2 -\-Pi < 1 ^ I 5 

4. 

h{i) 

Pi 

l-h{l)-Pi 

k <l <1% 

5. 

i-pp-Pi 

Pi 

Pf 

le < 1 


where 


l7{l)-=F- 


Cp - {y - cf) F {I - Pj - Pp) - CpF {I - Pj) 


Cf Cr 


□ 


Appendix E: Proof of Corollary [5] 

Proof. The proof of point 1 follows directly from the comparison of the closed-form dispatch 
solutions provided in Theorems [D and [2l For li>l 2 and 0 <l < min (p/ + ^i, ^ 5 ), dispatch rules 1-3 
in these two theorems apply and provide exactly the same solution {pw,Pf,Pi)- 

Statement 2 is a consequence of the fact that if pp >Pm, then I 2 = F~^ (^) Corollary 0]) , 

and thus, ‘"pfZ > implies that li>l2. Furthermore, if > p^^, it holds that Pj + h <15, because 

Cp — Cp Cp 

substituting pj + li into (fTBl) gives (u — cf) F {Ij —pp) -|- cfF{li) = Cp '^11 < Cp, since cf > Cp. 
From this point on, the proof of statement 2 proceeds as that of point 1. 

Claim 3 is based on Corollary 01 Indeed, the dispatch rule provided by this corollary is the same 
as that given by Theorem 0] when I ^ {h,Pw +Pi+Pp)- Furthermore, note that Iq <Pj+Pp +Pwj 
since substituting l = Pj+Pp+Pw into dUI) yields {v - cf) F ifpj+pp+pj^-pp-pj) + cfFifpj + 
Pp +Pw-Pi) = {y- cf) F (p^) -h c~pF{pp +Pw) = v>cp. 










46 


Morales and Pineda: On the Inefficiency of the Merit Order in Forward Markets with Uncertain Supply 
Article submitted to Operations Research^ manuscript no. (Please, provide the manuscript number!) 


Claim 4 stems from claim 3 by just noticing that =Pw +Pi when Cp = Cp and pp > py^. 

Statement 5 trivially follows from the fact that the dispatch solution provided by Theorem [2] 
is one of the possible solutions that the stochastic dispatch rule admits when cp = = Cp (see 

Corollary [3]) . 

Finally, the last claim relies on the fact that if Cp = Cp, maximum cost-efficiency is always 
achieved by breaking the merit order (see Corollary [T]) except for l>li+Pj when l 3 >Pp, and I > 
Pp +Pj otherwise. Recall that the price-consistent conventional two-stage market always prompts 
dispatch solutions that respect the merit order. □ 


Appendix F: Proof of Theorem [3] 

Proof. We divide the forward dispatch solutions that are feasible for model (12X1) into two groups. 
The first group includes those solutions in which the inflexible generation is dispatched below 
its capacity {pi <Pi). The solutions of the second group are characterized by the fact that the 
inflexible power generation technology is dispatched to its maximum capacity {pi = Pi). 

If Pj < pj and provided that Cj <Cp, the dispatch of the flexible generations^ has to be equal to 0 
to comply with the optimality condition of the lower-level problem fl21hh - (l21kh . Or, in other words, 
the merit-order forward dispatch imposed through (I21h() - (l21k(l implies that the more expensive 
flexible generation is only dispatched if the cheaper inflexible generation has reached its maximum 
capacity. For pp = 0,pj = l—pw and using the expected balancing cost of Proposition [21 the bilevel 
problem (I21h reduces to the following single-level optimization problem: 


Minimize 

pw>o 


s.t. pw < I 


rPw-PF rPw 

Ci{l—pw) + v / F{uj)duj + cf / F{oj)duj 

Jo JPw ~pf 


(48a) 

(48b) 


Based on the KKT conditions of this convex optimization problem, the dispatch rule and system 
marginal expected cost in this case is given by: 


Pw 

Pi 

Pf 

Marginal expected cost 

applies if 

1 

0 

0 

{v-c+)F{l-pp)^c+F{l) 

0<l<l2 

12 

I-I 2 

0 

Cl 

h <l <l2 +Pj 


Finally, the total expected cost corresponding to solutions with pj < Pj, which we denote by 
is computed as: 



where ^4(5) = {v — cf)F {s —pp) H- cfF {s). 


A{s)ds + Cj max(/ — I 2 ,0 ) 


(49) 
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On the other hand, the second group of feasible solutions are characterized by the fact that the 
inflexible generation technology is dispatched at full capacity, i.e., Pi =Pi- The dispatch rule of the 
flexible and stochastic generation can be thus derived from Proposition [2] with p = I —pj, that is. 


pw pi pF 

Marginal cost 

applies if 

l-Pi Pi 0 

h pj l — h—pj 

l-Pi-Pp Pi Pf 

{v-c+)F{l-pp)+c+F{l) 

cf + {v - Cp) F {1 -pp) - Cp (1 - F{1)) 

{v-Cp) F{l-pp) + CpF{l) 

Pi<l<h +Pi 

h +Pi <l<h +Pi +Pf 

h +Pi +Pf < 1 


Consequently, the total expected cost in this second case, which we denote by is given by: 


z^{l) = CiPj + / ^(s)ds+ / B{s)ds+ / C{s)ds (50) 

•y 0 J min{l—pj fi) J min{l—pj ,li-\-pp) 

where B{s) = Cp + {v — cj) F {s —pp) — Cp {1 — F {s)) and (^(s) = (a — Cp) F {s —pp) + CpF (s). 

Note that if l<Pj, the optimal forward dispatch must necessarily belong to the first group of 
feasible solutions and the total expected cost is, therefore, equal to z^{l). Likewise, if / > ^2 +P/, 
the optimal dispatch requires pi =Pi and hence, the total expected cost is given by z'^{l). For 
load levels pj < I < I 2 +Pj we must, however, compare z^{l) against z'^{l) to determine whether it 
is optimal to dispatch the inflexible generation at its maximum capacity or not. To conduct this 
comparison, we distinguish between two cases, namely, I 2 < h and li <l 2 - 
If I 2 <li, and given that Pj <l <l 2 +Pi, then I —pi <li and therefore: 


fmin(l-pj ,li+pp) 


l- 


pl-pi 

z^{l) = cipj+ / A{s)ds 

Jo 

Without any further assumption, we can rewrite z^{l) as: 


^min(Z,Z2) 


A{s)ds 


Pi-pj 

= C/max(/— ^2,0) + / ^(s)cis + 

Jo 

Therefore, for Pj < I < I 2 +Pi and I 2 we have: 

/ min(i,i 2 ) 

A{s)ds 

-Vt 


-PI 


(51) 


(52) 


(53) 


L-pi 


Next we evaluate the function z‘^{l) — z^{l) at the extremes of the interval Pj < I < I 2 + Pi and 
also compute its derivative with respect to I and obtain the following results: 




C 1 I 2 - /o^ A{s)ds > 0 

CiPi -ff Ms)ds>0 


if I 2 < Pi 
if Pi < I 2 


(54) 


z'^ih FPi) — z^{l2 APi) — 0 


( 55 ) 
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d {z^ {1) — {1)) )—Cl + A{1 — pfj <Q if l2<l 

m ~ \-A{l) + A{l-pj)<0 if iKh 


(56) 


where we have used that ^(s) < Cj for s < Z 2 by definition, and that ^(s) is an increasing function. 
Hence, we can conclude that z^{l) < z^{l) for Pj < I < I 2 +Pj and provide the following dispatch 
rule for l 2 <li- 


Pw 

Pi 

Pf 

applies if 

1 

0 

0 

0 <l<l2 

h 

l-h 

0 

h < 1 < h +Pj 

i-Pi 

Pi 

0 

b +Pi <l <li +Pi 

h 

Pi 

l-Pj-h 

b +Pi <l<h +Pi +Pf 

i-Pi-Pp 

Pi 

Pf 

h +Pi +Pf < ^ 


Let us now analyze the case h < I 2 for Pj <l <l 2 +Pi- First we rewrite z^{l) as: 

r 

z^{l) = CimeLx{l — l2,0) + / 

Jo 


+ / A{s)ds 

Jl-Pl 

And therefore: 


min{l—pj,li) pmin{l—pj ,li+pp) fl—Pl 

A(s)cis+ / A(s)ds+ / A(s)ds+ 

0 J min{l—pj fi) J min{l—pj ,li-\-pjp) 

(57) 




^m 

z^{l) — z^{l) = Cl {pj — m.ax{l — l2,0)) — / {A{s) — B{s))ds— 

J min{l—pj,li) 

~ /‘min(/,/ 2 ) 

{A{s) — C{s))ds — / A{s)ds 

) Ji—Pt 


(58) 


' niin(l—pj ,li+pp) 


By evaluating z^{l) — z^{l) at pj, I 2 +Pi and li +Pi we obtain: 

_ ^ _ J C 1 I 2 - A(s)ds >0 if l 2 <Pi 


min{l 2 ,li+pp) 


z^{Pi) - z\Pi) = 

Z^il2+Pl)-Z\l2+Pi) = - 

Jh 

z^{li+Pi)-z\li+pj) = 


CiPj - Jf A{s)ds>0 if Pi<l 2 


(59) 


{A{s) — B{s))ds—j' {A{s) — C{s))ds <0 (60) 

JTain(l2,h+PF) 


CiPi - A{s)ds > 0 


(61) 


if +P/ < I 2 

Ciih-h)-A{s)ds>0 if l 2 <h+Pi 
where we have used that A(s) > B{s),'\/s > li, that A(s) > C(s),Vs, and that A(s) < C/,Vs < l 2 - 
After checking that < q for < / < /a + P/, we can conclude that there must exist 

at least one OR THERE MUST EXIST A UNIQUE?? such that h+Pi <l 8 <l 2 + Pi and 
^^(/s) — z^{l8) = 0. Therefore, the dispatch rule if li < I 2 is: 
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Pw 

PI 

Pf 

applies if 

1 

0 

0 

0 <l < min(i2, h) 

h 

l-h 

0 

CO 

VI 

VI 

(N 

h 

Pi 

l — Pj — h 

ls<l< h-\-Pi -i-PF 

l-Pi-P f 

Pi 

Pf 

maxfisjh +Pi -\-Pf) < 1 


where 


h ■= X : li + pj < X <l 2 + Pj and C/ {pj — max(x — hjO)) — 


r‘mm{x—pj ,li-\-pp) 


f min{x—pj,li) 


(^(s) — B{s))ds— 


px-pi 


' rmn{x—pj,lx+pp) 


(^(s) — C{s))ds — 


^min(a:;,Z 2 ) 


A{s)ds = 0 


(62) 


-PI 


Consequently, the optimal dispatch rule prompted by the bilevel linear program f|2ip can be 
formulated as follows: 


Rule # 

Pw 

Pi 

Pf 

applies if 

1. 

1 

0 

0 


0 < / < /2 

2. 

h 

l-h 

0 


h ^ 1 ^ h -\-Pi 

3. 

i-Pi 

Pi 

0 

h < li 

h -\-Pi <l<li -hPi 

4. 

h 

Pi 

l-Pj-h 


h -\-Pj <l<li -\-Pj -hPF 

5. 

i-Pl-PF 

Pi 

Pf 


h -\-Pi +Pf < 1 

6. 

1 

0 

0 


0 <l < min(Z2,4) 

7. 

h 

l-h 

0 

h < h 

/2 < / < ^8 

8. 

h 

Pi 

l-Pj -h 


h<l<h -\-Pi-\-Pf 

9. 

I-Pi-Pf 

Pi 

Pf 


ma,x(l8,li -tpj +Pf) < ^ 


□ 


Appendix G: Proof of Corollary [6] 

Proof. Statement 1 is trivially inferred by comparing the tables provided in Theorems [T] and [3J 
rules 1-5 in both theorems are identical. These rules apply for li >12- 

Statement 2 follows from the fact that any stochastic dispatch solution that satisfies the merit 
order complies, by definition, with the following two conditions simultaneously: i) it is an optimal 
solution to the lower-level problem (I21hl) - (|21kp . because it respects the merit order and ii) it 
minimizes the expected system operating cost (|21al) . because it is a solution given by the stochastic 
dispatch rule. 
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Statement 3 simply highlights a particular case that is already covered by statement 2, since any 
dispatch solution for which pj =Pj preserves the merit order (recall that Ci <cp). 

Finally, statement 4 can be inferred from Corollary [3] and by noticing that, if Ci? = = c^, then 

^8 = ^2 + Pj- n 
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